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Recently a high precision spectroscopic investigation of the i5F^E^-X^E^ system of molecular 
hydrogen was reported yielding information on QED and relativistic effects in a sequence of rota- 
tional quantum states in the X^'T,'^ ground state of the H[2 molecule [E.J. Salumbides et al., Phys. 
Rev. Lett. 107, 043005 (2011)]. The present paper presents a more detailed description of the 
methods and results. Furthermore, the paper serves as a stepping stone towards a continuation of 
the previous study by extending the known level structure of the EF^'E'^ state to highly excited 
rovibrational levels through Doppler-free two photon spectroscopy. Based on combination differ- 
ences between vibrational levels in the ground state, and between three rotational branches (O, Q 
and 5* branches) assignments of excited EF^Y.'^ levels, involving high vibrational and rotational 
quantum numbers, can be unambiguously made. For the higher EF^T,'^ levels, where no combina- 
tion differences are available, calculations were performed using the multi-channel quantum defect 
method, for a broad class of vibrational and rotational levels up to J = 19. These predictions were 
used for assigning high- J _EF-levels and are found to be accurate within 5 cm~^. 



I. INTRODUCTION 

The hydrogen molecule and its deuterated isotopomers 
is the benchmark system for testing quantum ab initio 
calculations of molecular structure at ever increasing ac- 
curacy. While the first calculation of the dissociation en- 
ergy of H2 by Heitler and London Q , as the first applica- 
tion of quantum mechanics in molecular physics, was off 
by some 30% from the experimental value, the accuracy 
has improved by many orders of magnitude since then. 
The calculations of Wolniewicz in the 1990s obtained 
an accuracy of 10^^ for the binding energies of rovibra- 
tional levels in the Hj ground state and long stood as a 
benchmark result |3j l2i] ■ In the past few years improved 
quantum ab initio calculations of the Born-Oppenheimer 
potential and of adiabatic and non-adiabatic corrections 
for rovibrational levels of the X^Y,'^ ground state of H2 
have become available Q • These were extended with de- 
tailed calculations of quantum electro-dynamical (QED) 
and relativistic effects to yield a theoretical value for 
the dissociation limit P, and for binding energies of all 
rovibrational levels in the X^Y,'^ ground state of the H2 
molecule [6] with an accuracy at the 10~® scale. 

These highly accurate level calculations were tested, in 
fact prior to the publication of the theoretical studies, by 
accurate measurements of the ionization potential (IP) 
of H2 and D2 A measurement of the IP of HD 
followed [sl so that there are now tests for the advanced 
QED theory of the hydrogen molecule for all three stable 
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isotopomers [10|. Previously, tests of QED in molecules 
had been restricted to one-electron systems like the HD+ 
molecular ion [ll| . Note that in the IP measurements the 
absolute binding energy of the lowest (J — 0) quantum 
state of parahydrogen is probed with respect to the en- 
ergy of the ion. The contribution of QED and relativistic 
effects (hereinafter jointly referred to as QED effects) to 
the binding energy of the lowest rotational level in the 
specific case of II2 is 0.7282 (10) cm^^ 6]; the accuracy 
of this calculation is implicitly tested in the IP measure- 
ments. 

The availability of highly accurate QED calculations 
formed the motivation for testing these phenomena in 
a sequence of rotational states in II2, thereby seeking 
to reach high rotational quantum numbers. A predic- 
tion was made that the QED contributions to the bind- 
ing energy should depend on v and J quantum numbers 
in an experimentally detectable amount (61]. As for a 
possible detection strategy, the level structure of the H2 
X^T,'^ ground state can be probed directly through its 
purely rotational spectrum [l2| and through its vibra- 
tional spectrum [l3|, but due to the weak quadrupole 
nature of the transitions such measurements have not 
been performed in molecular beams with Doppler-free 
techniques. Recent Doppler-limited investigations using 
cavity enhanced techniques reached a relative accuracy of 
10~^, the most accurate determination of the quadrupole 
transitions to date [l3, [Hi • These studies are limited to 
rotational levels J < 5 and pressure shifts need to be 
corrected for. 

In the present study we adopted the strategy of prob- 
ing the ground state level structure via electronic transi- 
tions. In view of its lifetime in excess of 100 ns [la| the 
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FIG. 1. Potential energy curves of the X^T,'^ ground state 
and the EF^T,'^ state of H2, as well as the the X^T,'^ ground 
state of the ion. The excitation scheme is indicated with 
arrows. The first four vibrations of the inner well are indicated 
with solid lines, while the first six vibrations of the outer well 
are indicated with dashed lines. The rotational structure of 
the inner well vibrations is shown in the inset where the J 
quantum number is indicated. 



EF^Y,'^ state is a logical target and sensitive, high pre- 
cision, Doppler-free spectroscopy of the EF^I^t-X^Y^t 
system has been amply demonstrated over the years [ITI- 
|21| . This scheme involving two-photon excitation from 
the X^T,'^ ground state is illustrated in Fig.[T]and allows 
us to determine accurate level energies, as well as spac- 
ings between rovibrational levels, in the ground state of 
the H2 molecule, up to high rotational quantum states. 
Initial results of this work have been published in a recent 
Letter 

Using the present two-photon scheme, a limiting factor 
for further tests of binding energies in the X-^T,'^ ground 
state is the unknown level structure of the EF^Y]^ state 
at high rotational quanta, J > 5. Transitions in the 
EF — X system can be assigned as long as the level ener- 
gies in the EF^Ti'^ state are known, either from theory or 
from experiment. Yu and Dressier [23j have determined 
level energies for EF{v, J) levels from coupled-channel 
calculations involving a Born-Oppenheimer potential for 
the EF-sta.te [13], adiabatic and non-adiabatic correc- 
tions as well as relativistic corrections ,2i|j. These 
level energy calculations are limited to J < 5 and dis- 
agree with present measurements by a few 0.1 cm~^ for 
the lowest vibrational levels, and increase for higher vi- 
brational levels. Up to £^^(15) (at 110 000 cm"! ) the 
deviations between observed and calculated levels remain 
below 2.5 cm^^. 

As part of the present study calculations of level ener- 
gies were performed based on the multichannel quantum 
defect (MQDT) method 27[ . Previously results were re- 
ported on EF{v, J) level energies for J < 5 and the cal- 
culations are now extended to J = 19. The results are 



used to assign the EF levels, in particular, the high- J 
levels in the various EF — X bands observed. 



On the experimental side a comprehensive high- 
resolution Fourier-transform spectroscopic study was 
performed in the visible and near-infrared range access- 
ing a broad range of rovibrational levels in the £'F^E+ 
state 28] up to u = 28. Of relevance for the present 
study are determinations of level energies up to J = 12 
for E{v = 0) and up to J = 5 for E{v ^ 1 - 3); it 
should be noted that the absolute accuracies of level en- 
ergies in Ref. [2^ depend on the determination of anchor 
lines in the EF state by the two-photon laser experi- 
ments [U, [2^ . Those levels predominantly localised in 
the inner well are labelled as E{v), while those localised 
in the outer well are labelled as F{v), where v enumerates 
the levels in each well, beginning at u = 0. Above the 
potential energy barrier between the two wells of the EF 
state such labelling loses meaning and levels are referred 
to as EF(v) where, for each J, v enumerates all levels in 
energy order, beginning at u = 0. 



Calculations of QED effects in the ground state can 
also be tested by measuring combination differences, i.e., 
rotational energy splittings of A J = 2 in the ^^5]+ 
ground state from pairs of 0{J + 2) and <5( J), and pairs 
of Q{J) and S{J — 2) branch transitions. Still, such a 
comparison requires an unambiguous assignment of two- 
photon transitions in the EF — X system involving states 
with high rotational quantum numbers. Therefore the 
present study targets two scientific issues: (i) to test 
QED calculations of = 0, J) ground state rovi- 

brational levels, and (ii) to provide unambiguous assign- 
ments for transitions in the EF — X system involving 
high rotational quantum numbers for v > Q. 



Interaction Chamber 




Phosphor 



MCP 




PDA or PDL Laser 
System 




FIG. 2. Schematic layout of the experimental arrangement 
used in the experiments. 
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II. EXPERIMENTAL DETAILS 

A schematic representation of the general experimental 
arrangement is shown in Fig. [2l The precision metrology 
experiments were performed with a narrow-bandwidth 
pulsed dye amplifier (PDA) laser system. It consists of 
a three-stage traveling-wave optical amplifier, pumped 
by a Q-switched Nd:YAG (yttrium aluminium garnet 
) laser, and is seeded by the output of a continuous 
wave (cw) ring dye laser. The system delivers nearly 
Fourier-transform-limited pulses of 50 mJ and 5 ns du- 
ration, at a bandwidth of ^ 100 MHz. The instru- 
ment is tunable in the range 610-660 nm, while running 
on [2- [2- [4- (dimethylamino)phenyl] ethenyl] -6-methyl-4H- 
pyran-4-ylidene]-propanedinitrile (DCM) dye, although 
continuous scanning can only be accomplished over 1 
cm~^ intervals. A detailed description of this PDA sys- 
tem and its application in molecular spectroscopy is given 
in Ref. (sO). A second laser used is a pulsed dye laser 
(Quanta Ray PDL-3), also pumped by a Q-switched 
Nd:YAG laser. This laser features a grating-based oscil- 
lator delivering a bandwidth of ^^0.06 cm~^ or ^2 GHz. 
The continuous tunability over a broad frequency span 
makes this system suited for recording survey spectra. 
The laser was run on DCM laser dye and was tuned across 
the entire bandwidth of the dye with the aim of recording 
all detectable EF — X transitions in this range. 

The pulsed output of both PDA and PDL lasers is up- 
converted in two stages to deliver pulsed radiation in the 
204-220 nm range. Frequency doubling is accomplished 
in a KDP (potassium dihydrogen phosphate) crystal, and 
subsequent mixing in a BBO (/3-barium-borate) crystal 
results in the third harmonic of the visible radiation. 
Typical output powers produced in the deep-UV are 0.4 
mJ/pulse. 

Frequency calibration of the lasers and determination 
of the transition frequencies is achieved by simultane- 
ous recording of reference spectra of molecular iodine. 
For the PDL system a sufficient number of lines in a 
linear (Doppler-broadened) I2 absorption spectrum are 
recorded in each scan to linearize the frequency scale and 
derive an absolute calibration sj]. For precision mea- 
surement with the PDA system saturated I2 absorption 
spectra are recorded using the cw seed-laser; the absolute 
frequencies of the hyperfine components are calibrated to 
1 MHz accuracy §3, [11] . A linearized frequency scale in 
the PDA-precision measurements is produced from the 
transmission fringes of an etalon (free spectral range = 
148.96 MHz), which is length stabilized by locking to a 
frequency stabilized HeNe laser. 

Two-photon Doppler-free measurements of the EF~X 
transitions are recorded using two counter-propagating 
deep-UV laser beams, both in the case of the preci- 
sion measurements and for the survey spectra with the 
PDL-based laser. The UV-beam is split in two and ar- 
ranged co-linearly as illustrated in Fig. [5] exact counter- 
propagation is achieved by aligning the beams as part 
of a Sagnac interferometer [34], to avoid Doppler-shifts. 



The high resolution PDA source is mildly focused with 
a single 1 m lens to produce sufficient intensity to excite 
the EF — X transitions, but at the same time avoiding 
ac-Stark effects (see Sec. llVl for further details). For the 
PDL system a dual lens setup consisting of two 25 cm 
lenses is used for increased intensity at the interaction 
point to excite weaker transitions. 

H^ ions produced via 2+1 resonantly enhanced multi- 
photon ionization (REMPI) are detected and recorded 
as the laser is tuned. For the precision measurements an 
auxiliary 355 nm laser is used to ionize in a 2-1-1' REMPI 
scheme, i.e. a two-color scheme. Pulse delay of this ion- 
ization laser by 30 ns helps to avoid ac-Stark effects. For 
the survey spectra obtained with the PDL system a one- 
color scheme is used. Ions are accelerated in a time-of- 
flight tube, 50 cm in length, towards a multichannel plate 
(MCP). The signal on the MCP is converted by a phos- 
phor screen and a photo-multiplier tube, gated in a box- 
car and digitally stored. The voltages on the extraction 
plates are pulsed to measure under field-free conditions, 
thereby avoiding dc-Stark effects. 

In order to produce high rotational angular momen- 
tum states (and vibrational excitation up to w = 3) of 
the H2 molecule a sequence of reactions is used, occur- 
ring in a beam of HBr molecules. A pulsed solenoid valve 
(General Valve, Series 9) in the source chamber releases a 
pulse of HBr molecules that enters through a (1.5 mm di- 
ameter) skimmer into a differentially pumped interaction 
chamber. In a first photolysis reaction HBr molecules are 
photolyzed by the same UV laser beam (in the range 204- 
220 nm), which is also used as the spectroscopy laser: 

HBr + hv^Y ^ H + Br (1) 

resulting in H atoms with a high kinetic energy of up to 
2 eV, which is then sufficient to overcome the reaction 
barrier 35] to undergo a secondary exothermic reaction: 

HBr + H ^ H2 + Br (2) 

The reaction dynamics of such processes have been inves- 
tigated in detail and it is well established that these re- 
sult in product internal state distributions with rotation- 
ally and vibrationally excited H2(ti, J) molecules [36l - [39l |. 
Heck and co-workers [40| used a similar reaction dynam- 
ics scheme, using DI as a precursor gas, with the explicit 
goal to determine transition frequencies involving high- J 
levels (J — 26) in D2; the latter study was not performed 
under Doppler-free conditions. 

As an example, in Fig. [3] recordings with the PDL sys- 
tem of the EF — X (0,0) (3(13) transition are displayed. 
One of the measurements is recorded with a single laser 
beam yielding a Doppler-broadened profile 0.9 cm~^ 
width) , while a second measurement was performed with 
two counter-propagating beams resulting in a Doppler- 
free spectrum; the recordings are area normalized to 
one. Both measurements were conducted under similar 
focussing conditions, at intensity levels where ac-Stark 
broadening is not significant for these linewidths. The 
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FIG. 3. (Color online) Measurement of the Q(13) EF'^Ej- 
X^E^ (0,0) two-photon transition with the PDL system un- 
der Doppler-broadened and Doppler-free conditions. For de- 
tails see text. 



width of the Doppler-free signal (^ 0.3 cm^^) is caused 
by the laser line- width, accounting for frequency tripling 
and two-photon excitation. Figure [3] shows the advan- 
tage of the Doppler-free geometry in the enhancement of 
signal strength. Due to the large coUisional energy of the 
H atoms (2 eV), H2 products are formed with a spread of 
translational energies resulting in the Doppler-broadened 
line shape. 



III. MQDT CALCULATIONS 

The assignments of the EF energy levels were made 
on the basis of first-principles non-adiabatic rovibronic 
MQDT calculations. These were carried out exactly as 
in a previous study pTj . by use of the same theoretical 
approach [4l|, and the same input parameters for the 
computations. 



Briefly, in Ref. 27[ bond length- and energy-dependent 
quantum defect matrices fj,{E,R) were extracted from 
the then best available ab initio clamped nuclei potential 
energy curves of Refs. [13, 111,13] ■ Three distinct sets of 
matrices, corresponding to ^Sg (3x3 matrix), ^Hg (2x2 
matrix), and ^Ag (1x1 matrix) molecular symmetry, 
respectively, were used, each composed of singly excited 
lagciX and doubly excited lcr„e^'A Rydberg channels. 
IcTg and 1(7,1 denote the ground state and first excited 
state core orbital, and the partial wave indices £X 
and £'X' {£X, £'X' < 2) for the Rydberg electron are chosen 
such as to yield the required molecular symmetry. Each 
quantum defect matrix element is a smooth function of 
bond length and energy. The matrix elements were ad- 
justed to the quantum-chemical potential energy curves 
in such a way that when re- injected into a clamped-nuclei 
MQDT calculation, they reproduced the potential energy 
curves as best as possible. 



Rovibrational motion and ro-vibronic interactions were 
introduced by means of the frame transformation de- 
scribed in Ref. This transformation converts the 
clamped-nuclei quantum defect matrices into much larger 
rovibronic quantum defect matrices which typically are 
of dimension « 400 x 400 and which accoimt for the elec- 
tronic as well as for the rovibrational degrees of free- 
dom, and which, in particular, also include non-adiabatic 
electron-core interactions. Quantum defect techniques 
described in Refs. [13, Si], [l^] yield the desired non- 
adiabatic level energies. This purely ab initio approach 
reproduced more than 270 excited-state singlet gerade 
levels of H2 in the range < J < 5 with an overall rms 
error of about 6 cm~^ [23|, just slightly larger than the 
rms error of about 4 cm"-'^ obtained simultaneously by a 
coupled differential equations method for the same set of 
levels [13. 

This accuracy, which at the time corresponded to the 
state of the art of theory, is sufficient for the present pur- 
pose of assigning the new levels with J > 5, and therefore 
we have used the old approach without change. Improve- 
ments are possible but deferred at this time. They would 
include: (i) use of the improved first-principles potential 
energy curves available today; for instance, the poten- 
tial energy curve for the i?i^^S+ state reported in Ref. 
[25! and not accessible to the authors of Ref. [l^l is low- 
ered with respect to the curve used in [l^l by 1 up to 
15 cm^^, depending on the i?— value; (ii) an improved 
fit of the quantum-chemical data by the clamped-nuclei 
quantum defect matrices; indeed, the matrices derived 
in Ref. [131 reproduced the ab initio potential energy 
curves only to within about 4 to 8 cm"'^; and (iii) inclu- 
sion of a larger number of channels beyond t^£' < 2 in 
the clamped-nuclei MQDT treatment. 

Figure m shows all levels in the range 99 000 to 110 000 
cm~^ for < J < 19. The corresponding numerical 
information, along with existing high-precision experi- 
mental data are presented in Table I of the supplemental 
material (45|. The figure displays the MQDT level en- 



ergies and, as detailed in the caption, an indication as 
to whether the levels were seen in the present work, the 
work of Bailly et al. [1^, or both. Shading of the en- 
ergy levels in the figure indicates the character of the 
MQDT wave function for each level, either s-character 
(blue shading online) or combined doubly-excited and d- 
character (light pink shading online). The potential en- 
ergy barrier between the inner and outer wells, including 
the centrifugal J(J -I- 1) term, is shown by the essen- 
tially straight thick line (grey online) rising up from near 
104500 cm~^ at J = 0. In the region near and below this 
potential energy barrier those levels with predominantly 
s-character nicely line up to form the Eiv) inner well 
vibrational states. These levels are connected by wide 
shaded lines (light blue online) and are labelled i?(0), 
E{1), and E{2). Well below the potential energy barrier 
levels with predominantly doubly excited character line 
up to form the F(v) outer well vibrational states, partic- 
ularly at higher J values. These are labelled F{0 — 7). In 
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FIG. 4. (Colour online) Calculated energies and electronic character of the rovibronic energy levels of the EF^'E'^ of H2 shown 
plotted vs. J. The type of square used to show a level indicates the level's status: Levels observed by Bailly et aZ. [2^ are 
shown by squares with color-indented corners (red online), those observed in the present work are shown by double squares 
with indented sides (black online), while unobserved levels calculated in the present work are shown by simple squares (blue 
online). As can be seen some levels were observed in both Ref. [2^ and the current work. The calculated electronic character of 
each level is indicated by the size of the shaded square areas inside each outer square: Heavy shading (blue online) indicates the 
percentage of (lag) esa Rydberg s channel, while light shading (pink online) indicates the amount of doubly excited (Ictu) 
epX (0 < A < 2) channels as well as of singly excited (IcTg) edX channels. A fully filled square corresponds to 100% character. 
To guide the eye, the rotational progressions -E(O), E{1) and E{2) associated with the inner potential well of the EE state are 
connected by shaded lines. Heavy shading (grey online) indicates the potential energy barrier and the EF{v = 13) level which 
is above the barrier. Note that on the scale of the figure the largest observed-MQDT deviations are less than one third of the 
thickness of the lines used to draw the axes. See text for more details. 



the region well above the potential energy barrier the lev- 
els have predominantly doubly excited character and are 
now easily identifiable as vibrational states of the com- 
bined EF well. One of these, EF (13), is labelled and 
connected by a line (grey online). Beginning at J = in 
the region 105 500-106000 cm~^, just above the potential 
energy barrier, a renewed concentration of s-character oc- 
curs. For each J this s-character is shared between two 
adjacent levels for low values of J. The lower of each of 
these pairs was identified as belonging to an "i?(3)" level 
in Ref. [28| , but since the s-character is well less than 50% 
in each case the physical assignment no longer holds, al- 
though we retain the E(3) in the nomenclature for the 



purpose of comparison with literature values. This con- 
centration of s-character above the barrier is reminiscent 
of a shape resonance embedded in a dissociation contin- 
uum, except that in the present example the continuum 
is replaced by the relatively dense discrete manifold of 
levels associated with the outer F state well. Another 
recurrence can be faintly seen at low J around 107500- 
108000 cm-i in Fig. [H 
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FIG. 5. The EF^T.+-X^T.+ (0,0) two photon Q(15) transition 
recorded with the high-resolution PDA system. The fringes 
of a frequency-stabihzed etalon with FSR = 148.96 MHz is 
used to linearise the scan while the known I2 hyperfine com- 
ponent, marked with a o, provides an absolute calibration. 
The transition frequency axis (lower axis) is exactly six-fold 
the fundamental frequency axis. 



IV. PRECISION METROLOGY 

Precision metrology results on 21 transitions in 
the EF^Y.+-X^Y.+ {Q,0) band were obtained via 2-^1' 
REMPI from measurements on a line-by-line basis with 
the high resolution PDA system. As an example a record- 
ing of the (5(15) two-photon line is shown in Fig. [51 where 
simultaneous recordings of markers of the frequency- 
stabilized etalon and the I2 saturated absorption spec- 
trum provide the calibration. The spectral intensity of 
the light source is sufficient to detect S{J) and 0{J) 
branch transitions, known to be considerably weaker 
than (5(J)-branch transitions [4^. The A J = 2 ground 
state rotational level spacings were determined by taking 
combination differences between Q{J) and S{J — 2), and 
Q{J) and 0{J + 2) transitions. 

The identification of the two-photon lines in the 
EF — X (0,0) band derives from two consistent meth- 
ods: (i) the derived AJ = 2 combination differences in 
X^S+, V = match the calculations of Komasa et al. Q; 
(ii) the observed transition frequencies match differences 
between the highly accurate (experimentally determined) 
i?F-level energies by Bailly et al. [1^ and the X-ground 
state level energies of Ref. [q. The EF^Y.+ ,v = 0, J = 13 
level of Ref. [2^ is found to be inconsistent with observed 
combination differences, indicating a mis-assignment in 
Ref. ^2811. The J = 13, J = 15, and J = 16 levels were 
unambiguously assigned and measured at high precision 
with the PDA system. In view of the limited tuning 
range of the PDA system the J = 14 and J = 17 levels 
were not found. A listing of transition frequencies, cor- 
rected for the ac-Stark effect, is given in Table II in the 
Supplemental Material 45]. 



Systematic investigations were pursued into the uncer- 
tainties of the transition frequencies measured with the 
PDA-system. These are summarised in Table U and the 
most significant contributions are the chirp-induced fre- 
quency shifts in the PDA, the ac-Stark or power-induced 
effects, statistical fitting, absolute calibration of the fre- 
quency scale, and Doppler effects. 

A major contribution to the measurement uncertainty 
derives from the frequency chirp due to time-dependent 
gain in the dye amplifiers. This phenomenon has been 
extensively characterized for the presently used PDA sys- 
tem [i^l- Based on those investigations we estimate an 
upper limit to a possible chirp-induced frequency shift. 
From measurements on the Lyman and Werner bands in 
H2 it was found that at the edge of dye tuning curves 
this systematic effect can be pronounced [l^. Hence 
we include a conservative limit of 150 MHz due to the 
chirp in the error budget. A measurement of the Q(5) 
EF'^Yi^-X'^Yi^ (0,0) transition was found to agree within 
10 MHz with a previously reported measurement taken 
at higher accuracy, using a laser system with Fourier- 
transform laser pulses at longer duration [29] , indicating 
that the estimate of the chirp is conservative. 

The ac-Stark effect was found to broaden, shift, and 
cause an asymmetry in the measured line shapes of the 
EF — X two-photon transitions. This phenomenon was 
assessed by measuring each transition in the EF^Y,'^- 
Ar^i;+ (0,0) band as a function of intensity of the spec- 
troscopy laser in the deep UV; as an example results for 
the frequency shift are shown in Fig.[S]for the (5(13) tran- 
sition. An extrapolation to zero input power through a 
weighted linear fit yields the un-shifted transition fre- 
quency. The j/-axis error bars are the statistical uncer- 
tainty in each measurement point (averaged over three 
recordings) and the a:-axis error bars represent the error 
in the average pulse energy measurement estimated at 




Relative laser intensity (arb. units) 

FIG. 6. The measured g(13) EF^T.+-X^T,+ (0,0) two photon 
transition frequency at several different input powers. An 
extrapolation of a linear fit yields the un-shifted transition 
frequency. 
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10 % of the measured value. The power measurement 
was done with a UV sensitive photodiode and provides a 
relative power scale. The uncertainty in the zero power 
extrapolation fell between 0.0005 and 0.0015 cm~^. We 
take the upper limit, corresponding to 45 MHz, as an 
estimate of the uncertainty in the ac-Stark shift. 

The precision measurements are performed at low laser 
power (200-400 /iJ/pulse) inducing the two-photon exci- 
tation, and ionization by a 355 nm pulsed laser, delayed 
by 30 ns. In the case of the Q(6) transition the ac-Stark 
effect was measured on an absolute power scale. It was 
found that for the input powers and focussing conditions 
of the PDA system (200-400 /iJ/pulse focussed with a 1 
m lens), the ac-Stark shift amounts to approximately 150 
MHz. The ac-Stark shift coefficient for the Q{6) transi- 
tion results in a value of 1 MHz/ (MW/cm-^) which a gree s 
with the more accurate result of Hannemann et al. [2l| . 

The excitation is performed under field-free conditions 
to prevent dc-Stark effects shifting the transition fre- 
quency. Pulsed voltages of 1465 V and 2000 V are applied 
to the extractor and repeller plates respectively in order 
to collect the ions produced from the 2+1 REMPI pro- 
cess. Nevertheless, stray fields may be present that can 
cause a shift in the transition frequency and a conserva- 
tive estimate of 10 MHz is included in the error budget. 

The absolute calibration, derived from a Doppler-free 
I2 spectrum and a linearisation procedure based on the 
transmission markers of the stabilized etalon, is esti- 
mated at ^5 MHz in the two-photon transition frequency. 
Doppler shifts, possibly caused by deviations from a per- 
fect alignment of counter-propaga ting laser beams, are 
avoided in a Sagnac configuration [34|; residual angular 
mismatch amounts to some 1 MHz in the calibration un- 
certainty [2l[. The profiles of the observed lines in the 
Doppler-free geometry are essentially Lorentzian and the 
observed widths (on average 300 MHz ) are determined 
by the harmonically converted line width of the PDA- 
system. Line fitting errors were assessed by fitting with 
both a Voigt and a Lorentzian function. Deviations be- 
tween fits with the different functions were found to be 
at the 5 MHz level. The statistical uncertainty repre- 
sents differences on a daily basis and was found to be 
at the 30 MHz level due to the laser linewidth and the 
signal-to-noise ratio. 

In Table U a summary of the error budget is given re- 
sulting in a combined uncertainty of 0.005 cm~^ or 160 
MHz by adding individual contributions in quadrature 
sum. 



V. TEST OF QED IN THE H2 GROUND STATE 

From a combination of the precision EF — X transition 
frequencies and the known EF level energies [2^ ground 
state X^T,'^ rotational level energies can be derived. The 
accuracy of these results enables us to derive corrections 
to the ground state rotational levels that go beyond the 
nonrelativistic energy contributions. Experimentally it is 



TABLE I. Error budget for the precision metrology measure- 
ments obtained with the high-resolution PDA system. See 
text for details. 





Uncertainty MHz 


Uncertainty cm ^ 




1 




ac-Stark 


45 


0.0015 


Statistical 


30 


0.001 


dc-Stark 


10 


0.0003 


Line fitting 


5 


0.0002 


I2 calibration 


5 


0.0002 


Etalon non-linearity 


5 


0.0002 


Residual Doppler 


1 


0.00003 


Combined error 


160 


0.005 



not possible to disentangle the various contributions to 
the rotational excitation energies. However, starting with 
the experimental EF energies and subtracting the most 
accurate ah initio nonrelativistic energies of the ground 
state, comprising of the Born-Oppenheimer, adiabatic, 
and nonadiabatic contributions [J|, QED and relativis- 
tic corrections could be derived. This procedure yields 
the (differential) QED and relativistic corrections to the 
binding energies in the A"^E+ ground state for the rota- 
tional sequence J = 2 — 16 with respect to the J = 
corrections. 

Alternatively, ground state rotational energy splittings 
can be derived by using experimentally determined com- 
bination differences between EF~X Q{J) and S'(J — 2), 
and Q{J) and 0( J-l-2) branch transitions. By combining 
transitions with a common excited state level it is pos- 
sible to measure ground state rotational level spacings 
separated by AJ = 2. In similar fashion, by subtract- 
ing the BO, adiabatic, and nonadiabatic effects we 
derive differential J-dependent experimental values for 
QED and relativistic effects in the ground state rota- 
tional splittings. 

There are several issues in these analyses of the QED 
and relativistic effects that bear mentioning. Although, 
the two methods are equivalent in principle, the first 
method mentioned is dependent on the accuracy of the 
EF level energies determined in Ref. [2^, which are in 
turn referenced to the two anchor lines as measured in 
the EF — X system [2l|. The second method probing 
the ground state rotational energy splitting is indepen- 
dent of the accuracy of the EF level energies, although 
prior unambiguous assignments of the specific EF lev- 
els is a pre-requisite. The experimentally derived QED 
and relativistic corrections obtained from both methods 
are referenced to the X^S^, v = 0,J = level in the 
ground state. This level exhibits a combined QED and 
relativistic shift of -1-0.7282 (10) cm~^ towards the dis- 
sociation limit 6], and the entire rotational manifold is 
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FIG. 7. A plot of the QED effects (including relativistic ef- 
fects) in the X^Y:'^ ground state of H2 as a function of rota- 
tional quantum number J. Curves represent calculated dif- 
ferential QED effects (with respect to the v — 0, J = level) 
for the four lowest v = — 3 levels as calculated by Komasa et 
al. The data points with plotted uncertainties represent 
the experimental data as described in the text. 



shifted by this amount in the absolute sense. The ac- 
curacy of this theoretical QED and relativistic contribu- 
tions has been tested in an entirely different experiment 
sensitive to the absolute binding energy of the X^'S'^ 
V = 0, J = level Jj. The present work tests differential 
QED and relativistic effects in a sequence of rotational 
levels, up to J = 16 in H2. 

The experimentally derived QED and relativistic cor- 
rections are in excellent agreement with the recent most 
accurate calculations of Komasa et al. f§\ which includes 
the lowest-order relativistic corrections and QED correc- 
tions up to a^-order, where the accuracy of the calcu- 
lation is limited by the estimated higher-order a* QED 
contribution. 

The results of the experimentally determined QED and 
relativistic corrections in the II2 ground state are pre- 
sented in Fig. [7] (experimental data is given in Table 11 of 
Ref. [111). Also included are the theoretical predictions 
for QED and relativistic effects in rotational sequences 
for vibrations v — — 3. At present the predictions [6] 
and measurements [2I] agree within the experimental un- 
certainty, including the results of Wolniewicz |3| which 
cover the limited range of quantum states below J — 10. 
The EF state can also be used for testing theory for 
rotational sequences of vibrationally excited levels of the 
electronic ground state, provided that v > levels can be 
populated. A limiting factor for precision measurements 
towards further tests of QED and relativistic effects is 
the difficulty in the assignment of the transitions to the 
highly excited rotational quantum states (J > 10) of the 
EF state. The recording of survey spectra as presented 
in the next section is meant to address this specific prob- 
lem. 



VI. SURVEY SPECTRA AND FURTHER 
ASSIGNMENT OF BF-LEVELS 

Two-photon survey spectra of the EF — X system were 
recorded in the excitation range 90 000 - 97 000 cm~^ 
with the PDL system in a Doppler-free geometry with 
counter-propagating laser beams. Figure [8] shows many 
of the > 100 observed Q-branch transitions in the EF—X 
system. Lines are assigned with the convention followed 
by Bailly et al. [H] indicating the EFlv = 0), EF{v = 3), 
and EF{v = 6) vibrational levels of the inner well as 
E{{)), E{1), and E{2) respectively. The EF{v = 9) level 
is referred to as E{3) following Ref. f28l| although this as- 
signment is ambiguous if the wave-function composition 
is considered (see Sec. Illll) . The weak part of the spec- 
trum in the rectangle indicated in Fig. |S] is enlarged in 
Fig. ini demonstrating the large dynamic range of signal 
strengths probed in this investigation. 

Figure [9] shows a narrow range of the spectrum cover- 
ing some ^500 cm~^ where Q and S-branch transitions 
as well as lines connecting to the F outer well state are 
observed, ^-branch transitions are notoriously weak, but 
could be observed with the PDL system in the (0,1) and 
(0,2) bands, which have favourable Frank-Condon fac- 
tors [38] . However, these 5-transitions were limited to 
low, odd J only, as it is enhanced by the ortho-para ra- 
tio. 

All measured transition frequencies from the survey 
spectra are accurate to within 0.1 cm~^ and are available 
with assignments in an electronic document as Table 111 
of the Supplemental Material ^Sj ■ The two-photon tran- 
sitions from the survey spectra were assigned via three 
methods: (i) by measuring ground state energy splittings 
between O, Q, and S lines in a single band, as well as 
by ground state vibrational splittings between Q-lines; 
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FIG. 8. (Colour online) Doppler-free (2-1-1) REMPI spectrum 
of the H2 EF^T,+-X^T,j{v',v") bands between 90 000 and 
97000 cm^^. All assigned transitions belong to the Q-branch 
of their respective band. The small rectangle indicates the 
range shown in Fig. [9l 
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FIG. 9. Detail of Fig.[H]showing several transitions to both the 
inner, E, and outer, F, wells. 5* branch transitions (AJ = 
2) were observed for bands with favourable Franck-Condon 
overlap, in this case the E{0)-X{1) band. 



(ii) by making use of the highly accurate EF^Y,'^ level 
energies from Bailly et al. [2^ and the ground state level 
energies @; and (iii) by comparing to results from the 
present MQDT calculations. 

The E{0) and E{1) level energies were confirmed ex- 
perimentally by measuring ground state rovibronic en- 
ergy splittings, using both the PDL and PDA systems, 
and comparing to the calculations of Komasa et at 6] . 
For low J levels belonging to the E(2), E{5) and F(4) vi- 
brations the combination of the level energies measured 
by Bailly et al. [11] and the ground state calculations 
leads to an unambiguous assignment. 

A total of 30 transitions remained and were assigned 
by using the level energies from the MQDT calcula- 
tion and the ground state calculations. For each of 
these transitions the ground state rotational energies for 
ti = — 3 were added systematically to produce a level en- 
ergy which was then compared to the MQDT calculation. 
Levels matching within about 5 cm~^ were considered. 
Inspection of Table 1 of the Supplemental Material indi- 
cates that (with three exceptions labelled "tentative" in 
the table) the calculated wave function associated with 
each of these energies contains significant s and /or c? Ry- 
dberg channel inner-well character, thus making them 
accessible in vibronic transitions from the ground state. 
This feature provides further support for the correctness 
of the assignments. The majority of these transitions 
belonged to high J states with strong F character oc- 
curring near the crossings between E and F rotational 
levels as is indicated in Fig. |4l Furthermore high- J lev- 
els belonging to the E{2), E{3), and EF{10) levels were 
also observed. Finally, of these 30 transitions, four ad- 
ditional transitions correspond to upper-state levels that 
lie above the potential energy barrier between the inner 
and outer wells; these are indicated in Table HIl 



VII. EF^i:+ LEVEL ENERGIES 

Tables HI] and IIIII present the EE -level energies, for 
the inner and outer wells, respectively, determined from 
adding the rovibrational excitation energies ^ to the 
measured transition frequencies. Where possible, result- 
ing level energies (and their uncertainties) are averaged 
over various measurements of O, Q and S branches as 
well as over lines in vibrational bands. Included in Ta- 
ble |TT] are some highly accurate levels for w = and J < 5 
determined from a study with a more narrowband laser 
source [2l|, , as well as the accurate level energies de- 
termined by Bailly et al. (28j . 

When comparing a total of 38 measured £'i^^I]+ en- 
ergy levels obtained with the PDL system to those ob- 
tained by Bailly et al. [l^, a systematic deviation of 
-1-0.085 cm~^ between the data sets is found, which is 
attributed to an ac-Stark shift effect. In the metrology 
measurements with the PDA system it was found that for 
input energies of between 200 and 400 the ac-Stark 
shift amounts to 0.005 cm~^. The PDL-based survey 
measurements were performed at an intensity in the fo- 
cus increased by a factor of 16, which is consistent with 
an ac-Stark shift of 0.08 cm~^. The level energies pre- 
sented in Table |TT] recorded with the PDL system have 
been corrected for this ac-Stark shift. 

We find a difference of 2.95 cm~^ between our value 
for the E(l) J = 6 level energy and the value reported 
in Ref. The measured level energy derives from the 
E-X (1,1) Q(6) and E-X (1,0) Q(6) transitions, which 
in combination give an unambiguous identification. This 
strongly suggests a mis- assignment for the E{1) J = 6 
level in Ref. [Hi, in addition to the £'(0) J = 13 level 
which was also shown to be a misassignment 22]. Fur- 
thermore the assignment of the F{3) J=7 level differs by 
14.43 cm~^. This transition is based on the F — X (3,1) 
Q(7) transition which is assigned based on the MQDT 
calculations. Since the MQDT calculations are accurate 
to within ~5 cm~^ this is suggestive of a misassignment 
although experimental verification is needed to confirm 
this. In the present experiment this is not possible and 
hence we mark this assignment as tentative. 

The MQDT calculations of level energies were used 
for the identification of those levels which were neither 
present in the study of Bailly et al. [11] , nor observed in 
a combination difference. Figure [TUl shows the difference 
between experimental and MQDT-calculated level ener- 
gies. Overall these are < 5 cm~^. It had been found in 
the previous Ref. '^T] that while the MQDT calculations 
typically deviated by a few cm~^ from the observed level 
positions, these residuals turned out to be more or less 
constant along the rotational progression associated with 
a given vibronic level (e.g. E{v') or EF{v')), with abrupt 
changes occurring only near avoided crossings between 
inner- well and outer- well levels (cf. Table II of Ref. [27|). 
As Fig. [TU] shows, this feature is also borne out rather 
clearly in the new extended data set listed in the present 
Table [TTl For the -F(4) J = 6 — 8 levels, comparisons 
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with the current MQDT results and those in Ref. [28| 
show large discrepancies, in the order of 15 cm"'^. These 
levels were not observed in the present study using the 
PDA or PDL systems but from the evidence presented 
in Fig. [To] this is suggestive of a misassignment, although 
further experimental evidence would be necessary to con- 
firm this. 



VIII. CONCLUSION AND OUTLOOK 

A precision metrology study was performed on the 
i?F-'^E+-A^I]+ system of H2 via high resolution two- 
photon Doppler-free spectroscopy. The data were re- 
duced to provide a test of the differential, rotational 
quantum state dependent, QED and relativistic effects 
in the A"^E+ ground state of the H2 molecule which were 
calculated by Komasa et al. The agreement is excel- 
lent with the theoretical predictions currently more pre- 
cise than our experiment. In addition survey spectra were 
recorded in order to obtain an overview of the measur- 
able transitions in the £^i^^S+-A^E+ system, for higher 
V and J quantum numbers. This has served to extend 
the known level structure of the £^F^E+ state, thereby 
paving the way for future studies on QED and relativistic 
effects in vibrationally excited levels of H2 . QED calcu- 
lations show a net decrease for an increase in vibration, 
with effects of ^0.1 cm^^, well within the achievable ac- 
curacy of the PDA system. First-principles multi-channel 
quantum defect calculations were performed on EY{v, J) 
levels for high angular momentum states up to J = 19. 



The calculated level energies were found to deviate less 
than ~5 cm"-'^ from the experimentally obtained values 
and allowed the identification of otherwise unassignable 
transitions. 
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TABLE II. Table of level energies (in cm~^) for the four lowest vibrational levels belonging to the inner well and the first 
vibrational level above the barrier of the EF^T,f state derived from the measurements presented here and the ground state 
levels reported by Komasa et al. ^ . A represents the difference between the measured levels and the MQDT prediction. Results 
of Bailly et al. [23| are included for comparison. See the text for details on the levels grouped under tentative assignments. 
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101849.4044(2) 


4 


99791.32449 (10)'' 


3.57 


99791.32519(15) 102080.91 


y±'j) 


4.83 


102081.0311(2) 


5 


100098.26098 (Ibf 


3.55 


100098.26092(2) 


102367.17 


(\0) 


4.90 


102367.1451(2) 


6 


100461.196 (5)^= 


3.43 


100461.19733(2) 


102704.42 


y±'j) 


4.36 




7 


100877.369 (5)'= 


3.40 


100877.3708(2) 


103076.23 




1.39 




8 


101343.824 (5)'= 


3.40 


101343.82451(3) 


103525.37 




4.25 




9 


101857.174 (5)^= 


3.41 


101857.17482(2) 


103994.68 


(10) 


3.70 




10 


102414.045 (5)'= 


3.44 


102414.04588(5) 


104534.09 


(10) 


3.51 




11 


103010.497 (5)'= 


3.48 


103010.49993(3) 


105052.90 


(10) 


3.00 




12 


103641.544 (5)^= 


3.23 


103641.5428(8) 


105732.04 


(10) 


1.84 




13 


104307.453 (5)= 


3.60 




106232.97 


(10) 


2.13 




14 


105009.27 (10) 


1.75 




107004.11 


(10) 


2.59 




15 


105715.149 (5)^= 


3.62 




107494.19 


(10) 


1.23 




16 


106457.320 (5)^= 


3.85 




107917.60 


(10) 


-3.36 




17 


107186.40 (10) 


2.26 




108780.07 


(10) 


0.16 




19 


108777.75 (10) 


3.84 












J 


E{2) or 


EF{v 


' = 6) 


£■(3) or EF{v' = 9) 




Present Results 


A 


Bailly et al. [28] 


Present Results 


A 


Bailly et al. |28] 





103559.58 (10) 


0.97 


103559.59794(15) 


105384.90 


(10) 


0.42 


105384.9129(2) 


1 


103605.61 (10) 


0.46 


103605.6119(2) 


105415.28 


(10) 


0.23 


105415.2551(2) 


2 


103690.18 (10) 


-0.88 


103690.14695(14) 








105473.96704(15) 


3 


103995.28 (10) 


3.21 


103995.2119(2) 


105556.77 


(10) 


-0.85 


105556.8403(2) 


4 


104159.81 (10) 


3.74 


104159.80598(15) 








105657.7072(3) 


5 


104386.80 (10) 


2.74 


104386.8711(2) 


105770.05 


(10) 


-1.72 


105770.1314(3) 


6 


104650.43 (10) 


0.24 




105890.04 


(10) 


-1.80 




7 
9 


104908.69 (10) 
105826.53 (10) 


-3.10 
0.85 












11 


106894.85 (10) 


2.69 












15 


109293.90 (10) 


3.43 




107962.14 


(10) 


0.56 




16 


109607.40 (10) 


1.43 












17 


110526.86 (10) 


3.35 













J 



EF{v' = 10) 



Additional levels above the barrier 



Present Results A 

5 106374.02 (10) 3.52 

7 106721.14 (10) 1.55 

16 108937.93 (10) 0.73 



Bailly et al. [28] 
106374.1301(3) 



Present Results 
107537.79 (10) 
108386.73 (10) 
110202.77 (10) 
110806.31 (10) 



A 
8.35 
1.76 
4.64 
4.81 



Assignment 
EF{11) J=8 
EF{12) J=9 
EF{12) J=16 
EF{13) J=16 



^ Based on the measurements of Hannemann et al. 121 
Based on the measurements of Salumbides et al. [29 
Based on the measurements of Salumbides et al. 
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TABLE III. Table of level energies (in cm~^) for all measured levels belonging to the outer well derived from the measured 
transitions and the ground state calculations reported by Komasa et aZ. @]. A represents the difference between the measured 
levels and the MQDT prediction. The results of Bailly et al. [2^ are included for comparison. 



J 


F{3) or EF{v' = 5) 


J 


F(4) or EF{v' = 7) 


Present Results 


A Bailly et al. [28] 


Present Results 


A Bailly et al. [28| 


t 


103121.01 (10) 


0.62 


1 


103857.92 (10) 


-0.03 103857.8468(2) 




103420.55 (10) 


6.29 


2 


103903.07 (10) 


1.26 103902.9828(3) 




103672.48 (10) 


-0.91 


3 


103790.09 (10) 


-2.98 103789.9773(2) 








9 


104323.32 (10) 


-3.14 








14 


104973.02 (10) 


-2.29 


J 


F(5) or EF{v' = 8) 


J 




F(6) 



Present Results A Bailly et al. [28] Present Results A Bailly et al. [28] 
5 104971.99 (10) -2.05 104972.0087(3) 13 106727.84 (10) -1.91 

11 105523.80 (10) -3.19 14 106719.61 (10) -3.34 

12 105551.50 (10) -2.41 



* Tentative assignment. 



